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Abstract 

The response of inviscid incompressible unbounded fluid subject to a localized 
external perturbation is studied. The physically relevant hypotheses on the mode 
coupling mechanisms is justified by renormalization-group method. The scaling laws 
for the scalar and vector velocity potentials are derived. The spectrum of energy of 
perturbed fluid versus the distance apart from perturbation is computed. 

PACS number(s): 03.40 G, 47.55, 47.10. +g, 05.40. +j, 42.27. Gs 

1 Introduction 

The excitation of wave motion in a fluid by various perturbative factors is a long-standing 
problem in hydrodynamics. Particularly, the wave generation on the water surface by 
turbulent air flow over the surface and the rising of tsunamy-like waves excited by bottom 
earthquakes are of current interest of modern ocean studies. The physical mechanisms 
coupling the perturbation and water must be understood in order to describe the wave 
motions on the ocean surface. In spite of great success in the description of statistical 
properties of water waves which have been achieved recently in the framework of so-called 
Zakharov formalism (see, for example |I| and references therein) many phenomena on the 
ocean surface remain unclear. The great advances in study of surface phenomena were 
achieved in the modern critical phenomena theory recently 0. They are due to employ 
of quantum field theory methods and quantum field renormalization group method (RG) 
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0] in particular which are designed to describe systems with an infinite number of degrees 
of freedom. The results indicate convincingly that the analogous approach would also be 
of great advantage in the explanation of ocean surface phenomena. We intend to propose 
the statistical description of ocean surface phenomena within the framework of modern 
critical phenomena theory 0. In the present paper we consider a simple physical model 
coupling water to perturbation which can be adopted for the use of RG method. Our 
aim is to outline the basic properties of such a model for what purpose we now limit 
our consideration to the case of inviscid incompressible unbounded fiuid. The work is 
concerned with the large-distance and long-time asymptotic behavior of water response 
for an external localized perturbation risen by a pressure field pulse; it is imposed that 
the scale of perturbation Ip is much less than the main scale of the problem r 3> Ip. 
The duration r of energy input provided by the pressure pulse determines the scale of 
perturbed region 

Ip = CT, (1) 

where c is a speed of perturbation spread in fiuid. We are interested in the statistical 
properties of water response in the range r ^ Ip and t ^ t. The crucial importance 
for the what following is that the time derivative can be eliminated from hydrodynamic 
equations, 

div v(x, t) = 0, Ap(x,t) = -9iVj(x,t)ajVi(x,t), (2) 

where v(x, t) is the velocity of fiuid, p(x, t) is the field of pressure, and A is the Laplace 
differential operator. We have taken in (H) the constant fiuid density to be po = 1- 
Summation over successive indices will be implied, x Then from one obtains an 
expression for the pressure p(x, t) : 

JVp |x - y| 

where the integration is brought about over the perturbed region Vp. It is essential that 
(0) is invariant with respect to an arbitrary time-dependent velocity shift: 

Va(x,t) ^ v(x + s,t) -a(t), p(x,t) ^ p(x + s,t) (4) 

where a(t) is an arbitrary function of time decreasing at t — — oo, and s(t) = j'i^ a(t') dt'. 
This property expresses the Galilean invariancy of hydrodynamic equations. The conti- 
nuity equation in (^) shows that the velocity field v(x, t) can be presented as a sum of 
two terms 

v(x, t) = — grad 0(x, t) + rot A(x, t), (5) 
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where 0(x, t) and A(x, t) are the scalar and vector potentials consequently. The formu- 
lation of hydrodynamic equations in favor of the potentials if = {(f),Ai} elucidates the 
invariancy of the equations with respect to the shift of the vector potential Ai, 

A,(x)^ A,(x)-a,A(x) (6) 

in which A(x) is an arbitrary scalar function. This invariancy expresses the gauge symme- 
try of hydrodynamic equations. In the present paper we shall consider the consequences 
of this invariancy for the hydrodynamic equations (H) for the inviscid unbounded in- 
compressible fluid. Significant results on stationary spectra of the fully developed water 
response for the perturbative pulse may be derived phenomenologically in analogy with 
the fully developed turbulence theory and the statistical theory of waves based on the 
Zakharov's kinetic equations This approach is related to some conserved quantities, 
i.e., the wave action, energy, or momentum and to the idea of localization of pumping and 
dissipation ranges in separated parts of scale spectrum. By the way, a stationary spec- 
trum can have place for some scale interval which is transparent for the current of some 
quantity. For example, one can point out the well-known inertial range of Kolmogorov in 
the theory of fully developed turbulence which is transparent for the energy current from 
the large-scale region of energy pumping to the small-scale range of viscous dissipation. In 
the Kolmogorov's inertial range the energy distribution versus wave-numbers is described 
by the Five Thirds Kolmogorov's Law. In case of the fully developed turbulence the Law 
of Five Thirds describes the only possible spectrum since the three-dimensional isotropic 
movements of the inviscid incompressible fluid preserve the only integral quantity, i.e., the 
energy. However, in the problem discussed there are variety of spectra to be realized, since 
in the process of perturbation spread there are two more conserved integral quantities: 
the net momentum of fluid and the enstrophy. In case of an isotropic, 0(3)-symmetrical, 
perturbation the net momentum of fluid is equal to zero, however, even in this case there 
are two different spectra which are determined by the energy W— and the enstrophy 
currents. The velocity spectrum which is determined by the energy current from the 
pumping region kp = into the dissipation region A;^ = oo is the well-known spectrum of 
the fully developed turbulence theory, 

1^1/3^-1/3 ^ 

However, in the case considered fcp > fco — 7^ and kd ^ ko, so that the relevant 
inertial range lies apart from the scale spectrum ko > k ^ 1/r, and (|^ cannot have 
place in the problem considered. The enstrophy current (the squared averaged vorticity) 
determines the spectrum for the vector potential field A in the form 

S^/^k ^ A{k) (8) 
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which transparency interval is exactly the interval in question. In language of the critical 
phenomena theory the spectrum (^) determines the critical dimension A[A] of the field 
A, A[A] = 1. The use of critical dimension allows to compute the spectrum of any 
correlation function of the field A by simple dimensional counting. For example, for the 
pair correlation function in Fourier representation DA{k) = (A(k)A(— k)), one obtains 
an asymptotics: 

D^ik) ^ k^^''^\ A[DA] = 2A[A]-d, (9) 

where d is the dimension of space. For practical purposes, however, it is important to 
discuss not the correlation function (P) but a one-dimensional spectrum 

Da = ^X^'DA{k) (10) 

in which Sd = lix'^l'^ jVidj'X) is the square of surface of the unit sphere in dimensional 
space. For ([T0|) one obtains from (|^): 

Da ~ k. (11) 

Furthermore, the phenomenological considerations analogous to (|^) readily allows to de- 
termine the spectra of some quantities which can be measured in experiments, for example, 
for the energy as a function of distance from the perturbation point, E{r\ we obtain 

Er-^ ~ E{r). (12) 

From phenomenology, however, it is not possible to fix the spectra of all quantities si- 
multaneously. For example, one has nothing to say about the scalar potential spectrum 
(/)(fc). To determine it one can add some extra assumptions on the character of coupling 
mechanisms between different modes of fluid motions. Actually, the formulation of such a 
hypotheses is the crucial point of the problem of coupling water to perturbations. Clearly, 
these additional assumptions will sufficiently depend on the geometry of perturbation. As 
we have seen above, in some cases the perturbation pulse can rise a fiuid current, then 
the net momentum of fiuid P will be nontrivial. The new quantity conserved in the scale 
interval of the problem will leads to another possible spectrum. In the present paper we 
discuss the case of 0(3)— symmetrical perturbation; we consider the case with a current 
aligned in the fiuid in the forthcoming paper. Even for the simplest case of P = the 
formulation of physically relevant hypotheses on coupling water to perturbation is a non- 
trivial problem due to an infinite number of degrees of freedom. To solve it we apply the 
RG approach under the following reason: In the framework of RG method the physical 
degrees of freedom are to be replaced by the scaling degrees which are related to the phys- 
ical degrees through the RG transformations of fields and parameters of the theory. Since 



Critical Behavior. 



5 



the properties of scaling degrees of freedom possess a group structure (the renormahzation 
group) one investigates them much easier then those of the origin problem. The results 
obtained from RG-analysis are considered as somewhat statistical steady state limit of 
the physical system. The renormalized correlation functions are distinguished from their 
physical analogies only by normalization conditions, so that they can be also used for the 
analysis of asymptotic properties of the physical system. The plan of the paper follows: 
In the Section 2 we consider the direct consequences of the symmetry (^ for the hydrody- 
namic equations (^. As a result we derive an effective action functional which allows to 
understand the water response on the local perturbation as a critical phenomenon. The 
action functional presented has the most general form and then needs to be supplied by an 
assumption on the coupling mechanism. To formulate the physically relevant hypotheses 
we develop the RG-analysis in the Sections 3, 4, and 5 consequently. We show that in the 
symmetrical case considered (P = 0) the statistical properties of water response are man- 
aged by the dynamics of vector potential A. This conclusion instantly fixes the spectra 
of all quantities of the theory. In the Sections 6 and 7 we discuss the various composite 
operators of the theory which are responsible for the amendments to critical scaling and 
for the spectra of some quantities which can be measured experimentally. In particular, 
in the Section 7 we compute the one- dimensional energy spectrum E{r) as a function of 
distance from the point of perturbation. The results obtained for the spectra of the pair 
correlation function and the energy E meet the phenomenological relations ( [TT| ) and 
(|l^). We believe that the obtained results would give a key for better understanding of 
the origin of a numerous ocean phenomena such as Lengmuir circulations, spin-off eddies, 
and so on which theoretical justification have not given by this time. 



2 The Description of Model and the Effective Action 
Functional 

In this section we formulate a real time classical field theory of equations (^ in the 
Lagrange formalism and a model of nonlinear coupling mechanism between different wave 
modes. As a result we derive the effective action functional of abelian gauge-invariant 
field theory which describes the statistical properties of long-range water response on the 
external perturbation. We consider now the basic properties of the theory following from 
the symmetry (|^). 

(i) First, the eddy component of velocity field is expressed by a gauge invariant tensor. 



= (rot A)k = ^iAj - oijAi. 



(13) 
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(ii) Second, in accordance to the Noether's theorem the symmetry (j^) relates to a 
conserved current, vorticity, 

diJi = 0, dtJi = 0, (14) 

where Jj = djF^ - dkF^ , {i ^ k ^ j). 

(iii) Third, the classical equations (§) do not lead to a hamiltonian in the usual way: 
as a direct consequence of gauge invariance, the equations (|) do not depend on time 
derivative. It is therefore impossible to define the conjugated momenta, Nevertheless, 
(0) are just the field equations which can be derived from the classical lagrangian C{(p,p) : 



(dp)^ + {d<j))^ + ^F^ + p9iVj9jVi + JiAi 



(15) 



where we have introduced the tensor ViVj = {di(f)){dj(j)) + F^F^ — di(j)F^ — dj(j)F\ The 
statistical properties of mechanical system can be derived from the partition function of 
statistical mechanics Z = Tr{e~^) with somewhat classical dimensionless action S. In 
case of an infinite number of degrees of freedom one can write down the partition function 
in functional integral representation: 



Z{Ji) = J [d(p\[dp] exp —S{(p,p) + y ck dtJi(x, t)Ai(x, t) 
in which the euclidean static action S{Ai, 0) has the form 

S{y,,p) = \jdK ^dpf + {d<pf + ^Fj + p9iVj9jVi 



(16) 



(17) 



Here we note that the quadratic part of the action functional ([T7|) which is relevant to a 
free theory, i.e., with no coupling between water and perturbation is not symmetrical with 
respect to Galilean transformation (^ but is symmetrical with respect to the following 
transformations 

' A,(x)^A(x)-9.A(x), 
0(x) I — > 0(x)e^"«^W, 

where A(x) as usual is an arbitrary scalar function, and uq is a potential coupling constant. 
These transformations expresses the so-called U{1)— gauge symmetry (f/(l) is the group 
of multiplication by complex numbers). The last term in (|l^) does not meet the entire 
symmetry (p!8D but preserves the simple gauge symmetry (^. This symmetry breaking 
shows that the pressure field as it is included in the action functional contains somewhat 
redundant degrees of freedom. 

Instantly close to the region of perturbation the pressure pulse rises the wave motions 
with eigenmodes k > k^ ~ l//p- Due to strong nonlinearity of the interaction in the 
Navier-Stockes equation the eigenmodes of oscillations spread very fast from a band of 
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order over the whole spectrum, and various multipole oscillations of any type are 
arisen with time. One can say that after a short period of time motions of any modes are 
present in the water. Clearly, the long-distance fluid behavior will depend to some extent 
on the statistical properties of wave mode coupling. 
Assuming an effective action functional to be 

(i) local in space and time, i.e., it depends only on the fields A(a;) and (pix) and their 
partial derivatives (and not on products of fields and their derivatives at different points), 

(ii) invariant under space and time translations, i.e., space and time coordinates do 
not appear explicitly in the action, we now suppose the simplest model for the coupling 
mechanism by inclusion of the y^^-type interaction term into (p!7|): 



(19) 



with a wave modes coupling constant g. In accordance with the general critical phenomena 
approach we note that the accounting of highest oscillation harmonics, i.e., y?^, ip^ and 
so on cannot alter the large-distance asymptotic behavior of water response for localized 
pressure pulse if 7^ 0. 

Since the pressure pulse is localized in the scale / < /p, and the locality postulated for 
the action functional requires redundant degrees of freedom for the fluctuating pressure 
field p(x, t), we therefore can integrate it over in the partition function (pIBf ). The result of 
functional integration does not depend on p(x, t). Due to (^ this procedure is of perfect 
clarity and is reduced technically to elimination of the quadratic term proportional to 
[dpY from (|1^) and to replacement of the pressure field p{^,t) in the U{1)— breaking 
term pdiVjdjVi by (^. In particular, it leads to a new term in (|IU]) of the form 

^ f ^ » r 4-\^ f 4-\ f A ^iVj(y,t)^jVi(y,t) 

- / cix aiVj(x,t)ajVi(x,t) / dy — , _ , , (20) 

/ J JVp |x y| 

which relates the fluctuations of velocity fields risen by the perturbation pulse in Vp to 
the fluctuations apart from the perturbed region. We note that ( PUj) preserves the entire 
[/(I)— gauge symmetry (|T8|). 

Performing the integration over Vp and moving the derivatives onto the result of inte- 
gration, one can rearrange (pO]) to the form 

iy'rfxKij(x)Qij(y.,x,t), (21) 

where Qij is a quadratic form of potentials and Ai, and Kij is a source kernel. To explain 
the meaning of (^) we adduce some arguments from renormalization group analysis which 
we are going to apply to the theory discussed in the forthcoming sections. For the future 
purposes of RG analysis we need only insertions at zero momentum and we could, in 
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principle, restrict ourselves to constant sources Kij, since the renormalized action which 
we derive for the case in question is obtained by setting Kij{x.) as constant, [Q. 

Another important feature of (^I]) stems from the fact that the RG transformations 
generate all possible linearly independent quadratic terms in Qij, i.e., 0^(x, t), 0ylj(x, t), 
and AiAj{'x,t). However, some certain linear combinations of quadratic operators in Qij 
only are relevant to scaling degrees of freedom and possess, by the way, the definite 
physical meaning. We compute all such linear combinations in the Section 6. Nevertheless, 
one can see before calculations that the linear combinations of quadratic operators which 
are of importance for the RG analysis have to be f/(l)— gauge symmetrical. There is only 
one such a combination, i.e., 

m20'(x) (22) 

in which is somewhat mass parameter (the coefficient of the relevant RG-invariant 
operator). The use of Ward identities which express the f/(l)— gauge invariancy of the 
theory allows to show that all other combinations of quadratic operators are ultra-violet 
(UV) finite, i.e., the relevant correlation functions do not have UV divergences. Thus, with 
no loss of generality we can omit all combinations except for (|2^) from the consideration. 
It leads to the action functional which follows 



(23) 



The effective action functional (^) which is arranged to describe the asymptotic properties 
of water response is not renormalizable, and it has no solutions in the massless limit, when 

= 0,1^. Again, the reason is that it still has some redundant degrees of freedom, the 
gauge degrees, with unknown dynamics. 

To construct a renormalizable theory we are led to introduce an abelian gauge geo- 
metrical structure. By the way, 

(i) 0(x) and 0*(x) are vectors for f/(l) transformations, 

(ii) the derivative di is replaced by the covariant derivative V,: 

Vi = di + iuoAi, (24) 

where uq is the coupling constant of interaction between the scalar and rotational compo- 
nents of the velocity potential if (analogous to the electron charge e in electrodynamics). 

(iii) It follows that the curvature tensor is iu^Fij : 

iu^Fij = = iuo{diAj - djAi). 

(iv) Since the f/(l)-gauge group is abelian (A(x, t) is a translation invariant), one can 
write the parallel transporter U{C) along any continuous contour C which is an element 
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of U{1). In terms of a line integral: 



U{C) = exp 



-tuo (j> Ai{s) dsi 



(25) 



as a consequence of vorticity conservation. Thus, the rotational component of velocity 
potential just carries on the fluctuations of the scalar potential field 0(x). By the way, two 
solutions for different points 0(x, t) and 0(y, t) are related through the parallel transporter 
(p5|), where C is an integration path connecting the points x and y, 

The form of the action functional (|23|) which meets the geometrical structure (i)-(iv) 
follows 



(ix 



(5p)'|Vi^ 



+ 



(26) 



The problem of the gauge invariant theory is that the local gauge invariant action does 
not provide a dynamics to the gauge degrees of freedom. We take them into account, 
performing the standard procedure analogous to the Fadeev- Popov quantization [Q. Let 
us write the gauge field Ai in terms of a gauge field Bi projection of Ai on some gauge 
section, i.e., satisfying some gauge condition, and a gauge transformation: 



Ai = Bi + diA. 

We assume that this decomposition is unique. For A(x) one imposes: 

92A(x) + 9iS,(x) = /i(x). 



(27) 



(2^ 



in which h{x.) is a stochastic field for which a probability distribution is given. We do 
not include a term proportional to 0(x)e*^^'''' in the condition ([28| ) omitting the U{1) 
gauge degrees of freedom to simplify the model and to avoid the appearance of Fadeev- 
Popov ghost fields. Including the equation ( pSj ) in the functional integral for the partition 
function, one can see that since the result does not depend on the dynamics of A(x) and 
on the field h{x.) either, one can integrate over /i(x) with the gaussian measure: 



[dp{h)] = [dh] exp 



1 

2C 



-— / dch^i 



(29) 



in which ( is an arbitrary valued {( G [0, oo)) auxiliary gauge parameter of the theory. 
The resulting effective action of the model has the form: 



S{A,, 



dx 



r + ^4 + r^(5.A,)^+mV + ^#^ 



(30) 



In the forthcoming sections we show that the model (|3(]|) demonstrates the existence of a 
statistically steady state independently of the details of velocity evolution. The functional 
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(^) is analogous to a t/(l)— invariant action for a charged scalar field with a 101^^ self- 
interaction (so called Ahelian Higgs Model, [Q). This theory allows the multiplicative 
renormalization, [|10] and for some values of the parameters {mo, go, Uq} the corresponding 
physical system tends to a steady state in a large-distance limit and demonstrates a 
universal behavior. In language of the critical phenomena theory, there is a nontrivial 
IR-stable fixed point which determines the critical asymptotics of correlation functions of 
the fields 0(x) and A(x). The relevant critical index for A meets the phenomenological 
result d). 

We conclude this section by an explanation of physical meaning of solutions for differ- 
ent signatures of the mass parameter m^. In case of A = in (0), one has the standard 
model of a scalar unharmonic oscillator. This model is, may be, the most popular and 
well-investigated action of the modern theoretical physics. For m = the oscillator is 
subject to a phase transition. At the classical level in the symmetrical case (m^ > 0) 
the oscillator model describes the fluctuations having the trivial expectation value of the 
field, (0)o = (see Fig. l.a). If < 0, the system allows two possible expectation values 



for the field (0)o = ^\Jm? / g (see Fig. l.b). The latter situation is usually referred to as 
spontaneously broken symmetry. 

The physical consequences for the model considered can be readily understood. The 
signature of integral operator kernel i^jj(x, t) depends on the certain physical conditions 
and on the geometry of perturbation, and it is determined directly by the signature of 
operator diVjdjVi. This signature does not depend neither on the certain value of velocity 
nor its evolution, but it depends on the topological properties of fluid flow risen by the 
perturbation. 

The case of symmetrical perturbation, when the net momentum of fluid is equal to 
zero (for example, on the surface of a large scale eddy, see Fig. 2. a) can be described by 
the model (|30D with > 0. In the vicinity of saddle points, diVjdjVi < 0, i.e., when there 
is a net fluid current from the region of initial perturbation into outside (see Fig.2.b), and 
the net fluid momentum P ^ 0, one can use ( |30[ ) with < 0. 

In the what following sections of the present paper we shall consider the case of sym- 
metrical perturbation, > 0, (Fig. 2. a). The presence of the gauge field A affects 
sufficiently the behavior of the system due to the transporting role of the field A. We 
demonstrate below that the ordinary infrared (IR)-stable fixed point of RG transforma- 
tions which is responsible for self-similar behavior in the standard 0^-theory turns out to 
be unstable, however, the new fixed point acquires stability for the real value e,. = 1/2 of 
the parameter 2e = 4 — d. We shall derive formulae for universal scaling profiles to the 
first order in e. The profiles has a power-law behavior for the large distances r ^ Ip. 
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3 Infrared Singularities of Perturbation Theory Dia- 
grams 

In the present Section we develop the diagram technique relevant to the theory (|30|) 
and discuss the large-distance {k 0, in momentum representation) singularities of 
perturbation theory diagrams. 

The model (|30D can be considered in the ci-dimensional space x with UV-cut off A = 
ko — 1/lp- Each quantity in (|30D corresponds to one (momentum) canonical dimensionality 
df, which is completely determined by the space dimensionality d. 

In the critical phenomena theory one seeks the asymptotic for correlation functions in 
the region k, m <^ A for which one considers (70 — — const A^~'^ with const < 1. From 
now on we supply all the parameters in (|30|) by the lower index "0" to distinguish them 
from those in renormalized action forthcoming. In order to be specific we consider the 
pair correlation functions of the potentials D^{r) = (0(x)0(y)) and -D^(r) = (A(x)A(y)) 
in which r = |x — y|. 

In momentum representation these correlation functions are found from the Dyson 
equations {p is the external momentum) 

D^' =p' + ml- S^(p), Dj' =p' + ie- S^(p), (31) 

where ^^(p) and S^(j9) are the infinite sums of all 1-irreducible Feynman diagrams (see 
Fig. 3) whose vertices correspond to the multipliers go and Mq 




and whose lines correspond to the bare propagators 
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A(k) 



A{-k) 



where Pij{(,k)6ij + (C — l)kikj/k'^ is the gauge dependent projector, and e is a reg- 
ularization parameter for k = 0. We investigate the theory in = 4 — 2e dimensions 
considering e as a small parameter of a regular expansion which has 2er = 1 as an actual 
value. For < 2e < 1 the diagrams in Fig. 3, independent of p and mo, have the algebraic 
UV-divergent terms ~ ^2-^2ne (^-^j^gj^g ^ jg ^j^g order of perturbation theory) corresponding 
to the simple shift of rriQ which does not alter the signature of mg. If we consider the 
value of mQ to be known exactly, it is necessary to discard all such terms. As usual P], 
this is implemented by subtracting their values for p = tuq = from all graphs of Fig. 3. 
After these subtractions the integrals for < 2e < 1 become UV-convergent, the cut off 
A can be eliminated (taken as oo, i.e., the initial perturbation scale Ip is taken as zero), 
and the series Fig. 3 takes the form 



D- 



{p^ + ml) 



n,l=l 



(32) 



For p ~ mo -C A and e > the dimensionless parameter of the expansion gop' 



-2e 



{Ap 



-l\2e 



m 



is not small, and it is necessary to sum the series. This problem is solved 



by the RG method. 

In a clearer formulation this problem reduces to a determination of the asymptotic 
value of the propagator Dy, = DixViX^^o) for % — 5> (everything is fixed except for x). 
This procedure is nontrivial for e > due to presence in the c„^/ of poles in e and leads to 
the equations of RG which we shall consider in the next Section. 

Another problem, which is occurred in the region uiq <^ p, is connected with sin- 
gularities of the coefficient ; in (^) for vtlqIp and cannot be handled by RG. 
This problem originates from the finiteness of the physical value of e. After removing of 
UV divergences from diagrams of Fig. 3 there are still diagrams diverging for m — for 
any e > 0. This problem had been discussed in the critical phenomena theory where the 
method of short distance expansion (SDK) was employed, [Hj-pl. We shall apply SDK to 
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compute the leading amendments to critical scaling of ( pO] ) with > in the Section 6. 

Finally, we make a note on the particular features of perturbation series Fig. 3 for 
the gauge invariant theory. Some diagrams in Fig. 3 have ^"dependent poles which are 
unphysical, since they have been introduced to make the theory renormalizable. The 
renormalization constants of the gauge invariant theory as we shall define them later on 
are gauge independent, therefore we can fix the value of the gauge parameter ( in certain 
calculations. In particular, we use the Landau gauge = 0), so that the gauge field 
propagator is simply proportional to the transversal projector Pij. 



4 Renormalization- Group Equations. Scaling Degrees 
of Freedom 



Now we discuss the renormalization procedure for the model (|30| ) and produce the renor- 
malized action functional, then we derive the RG equations for renormalized correlation 
functions. Renormalizability of the theory (|30D (the Abelian Higgs Model) for any value 
of the gauge parameter ^ < oo, is proven (see for example ^^), and we do not discuss it 
in details. The UV-divergences (in our case the poles in e in diagrams) of the model con- 
sidered are removed by the multiplicative renormalization procedure. It amounts to the 
following: the initial action is referred to as nonrenormalized, its parameters and coupling 
constants are referred to as bare; these are considered as some functions (remaining to be 
determined) of new renormalized parameters and coupling constants. The renormalized 
action functional 



is a function of renormalized coupling constants and parameters: 

9q = l^'^^gZg, ul = fx^'u^Zu, ml = rn^Z, 



(33) 



mi 



2 _ 7 i2 a2 _ 7 . /12 



A' = ZaA'ri Co = ^cC/ 



(34) 



where all renormalization constants Za are the functions of four independent quantities 

Zi = Zfj,, Zu = Z(^ = Z = Z2 ^, ZmZ^ = Z3, Zi = ZgZ'^, (35) 

which can be calculated within the framework of diagram technique. We chose the simplest 
form of subtraction scheme, where the divergences are presented as the bare poles in e (so 
called "minimal subtraction scheme"); /i ~ l//p is the renormalization mass parameter, 
gXi'^ cind u are renormalized analogies of the bare parameters ^oCoj^o and Mq, Za = 
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Za{g,e,u,d) are the renormalization constants. Due to gauge invariance of the theory 
the terms breaking the gauge symmetry are not renormahzed, and they do not require 
counterterms, 0]. The renormahzed correlation functions Wr meet the relation 

W^{g, u, f^)Z^^Z^- = Wigo, (36) 

in which are UV-finite functions (they are finite in the limits e — > 0) for fixed param- 
eters a. 

The RG equations are written for the functions which differ from the initial W 
only by normalization and then can be used equally validly for critical scaling analysis. To 
derive these equations one notes that the requirement of eliminating singularities does not 
determine the functions Cq = eo(e,e), e = {g,u,m}, uniquely because of the value of /i 
is not fixed by any physical condition. Variation of fi for fixed values of bare parameters 
Co leads to variations of e and renormalization constants (|35|). Following the standard 
notation, we denote by the differential operator /U-D^ for fixed Cq. Applying it on both 
sides of (32) leads to the basic RG equation, 

[Df, + Pgdg + Mu - ImD^A = 0, (37) 
where we have used = xd^ for any parameters of the renormahzed theory; for any Zi 
7i = L'^lnZi, Pa = D^a, a = {g,u,}, i = {g,u,C,fn,(j), A}. (38) 



These identities determine the /^-functions of the theory considered. 



2e + 



(39) 



and the anomalous dimensionalities 7^. One calculates the renormalization constants Zi_4^ 
from the diagrams of perturbation theory (these calculations are completely analogous to 
the relevant computations in /^-theory of the critical phenomena theory, |^) and then, 
using (^), 7i— and /5a-functions. By the way, all 7— and /?— functions are constructed 
as series in g and u, and the functions 7j do not depend on e. Furthermore, the relations 
between renormalization constants (^) lead to analogous relations for 7^: 

7g = 74 - 271, Ju = -72- (40) 

We have computed the relevant renormalization constants up to the second order diagrams 
of perturbation theory. These computations are pretty standard, so that we just bring 
about the results for ■ji for the three dimensions, 

7i = V + «", 72 = 6n'2, ^^^^/2^Q^/2^ ^^ = Qg'2^1g' + 4u', (41) 

n .i 
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where g' = g/lQn'^ and u' = -u^/IGtt^. From (^) one obtains the exphcit expressions for 
the 7-functions of fields and the mass: 

70 = J^7' + «", 7A = 6M'^ 'ym = -9'-lu' + ^g" + 5u'^. (42) 
D -3d 

Substituting (^) into (^OD and (0), we obtains the expressions for the /3— functions: 

7 , . , 17 



/3„ = -u' 2e-6u , /3„ = -g 



2e+-g' + Au'+^g'^~2u'^ 



(43) 



Eight fixed points of the RG transformation are determined by the system of equations 
j3a{g'*,u'*) = 0. A fixed point is stable with respect to large-distance asymptotics if the 
matrix uJij = di(3j is a positively defined matrix at the fixed point. Fixed points and their 
stability regions with respect to the large-distance asymptotics are collected in the Tab. 
1. 

The fixed point 4 corresponds to the asymptotic behavior of the scalar 0^ model 
of critical phenomena, [0]. In the model discussed it is unstable due to scalar potential 
coupling to the vector potential field A. For the "physical" value of the parameter 2e = 
A — d the point A^ 3 is the only stable fixed point. The fixed point A^ 7 would be stable close 
to the four dimensions also. The large-distance asymptotic behavior of water response 
in three dimensions is then governed by the fixed point A^ 3. The inequality < e < 3 
determines the relevant basin of attraction in space dimensionality. 

In the framework of RG-approach the physical degrees of freedom are replaced by the 
scaling degrees including anomalies. In particular, the scaling functions are obtained in 
the form of the power series in g and u. Asymptotically, these coupling constants are 
replaced by their values in fixed points of RG-transformation. The critical indices in the 
4 — 2e expansion are obtained from the 7— functions (^) with replacement of g and u by 
g* and u* also. By the way, the properties of scaling degrees of freedom (see Tab.l) yield 
some qualitative conclusions on the physical properties of the model. 

From the data of Tab.l it follows that for A^ 3 the scaling degrees of freedom related 
to the 0^— coupling vanish {g* = 0) in three dimensions. Thus, all diagrams proportional 
to g are vanished in the symmetrical phase of (pO]), however, no one correlation function 
becomes trivial. Considering the theory in three dimensions, one can eliminate the scalar 
wave mode coupling term g(f)^ from the action (|30|) , since in case of > it does not 
contribute to the large- distance asymptotics of water response. One can say that three 
dimensional dynamics of an inviscid incompressible fluid involved in eddy motion has 
somewhat short of physical degrees of freedom to allow the coupling between the different 
scalar wave modes. A^ 7 in Tab.l gives us an evidence that the both coupling mechanism 
introduced in (pOF) are of equal importance in four dimensions. 
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Choosing the certain fixed point of RG-transformation {i.e., N 3), we neglect all cou- 
plings between various wave modes in benefit for the certain one which is responsible 
for vorticity conservation (|25|) in the large- distance limit. This conclusion expresses ex- 
actly that additional assertion which we have needed to complete the phenomenological 
description of the problem of coupling water to the 0(3)— symmetrical perturbation. 

5 Solution of RG-Equations. Critical Scaling and 
Asymptotics for Pair Correlation Functions 

In the present Section we derive the solutions of RG differential equations ( P7| ) for 
and Da- The use of standard dimensional counting supplied by (^) and (^21) leads to the 
following expressions for the one-dimensional spectra of the theory: 

D^{s) ^.^0 s'-'-''^^U{s,g,u',z), Da{s) ^.^o s'-'+'^-fA{s,g,u',z). (44) 

The functions and /a of dimensionless arguments: s = kip, z = m^r^/" , r]^ = (i — 2 + 7,^ 
= d-2 + e/3, riA = d - 2 + 'Ja = d - 2 + 2e, and 1/z/ = 2 + 7^ = 2 - 2e meet the 
RG-equations of the type (|37D which allow to find out their scaling asymptotics (s — 0). 

Here we note that for = 1/2 the asymptotics (0) for Da meets the phenomenological 
result ([Tl|). The spectrum for D^j) cannot be predicted from the bare dimensional counting 
and is justified within the framework of developed RG method. We have pictured these 
spectrum out in Fig. 4. 

Being the solutions of RG-equations, the functions and /a are to be the arbitrary 



functions of the first integrals of (|37|). The number of first integrals is one less than the 
number of arguments of / in (0), and they can be founded from the system of equations 

ds dg du ^^^^ 

s (3gig,u) /3ui9,u)' 
supplied by some normalization conditions for g and u. We use the standard one, 

g{s = l,g,u) = g, u{s = l,g,u) = u. (46) 

At the fixed point 3 from (|45|) and (|46| ) one obtains the asymptotic solutions for 

f<f>{.s,g,u^,z) and fAis,g,u^,m'^r^/''): 

f^{s,g,u)=[^] T^{l,u\z), (47) 

and 

f^{s,g,u) = Ta{1,u'\z), (4J 

the scaling functions JF^ and JF4 are not fixed by the RG-equations and calculated usually 
in the framework of diagram technique. 













r 
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6 Short Distance Expansion. On the Possible Cor- 
rections to Critical Spectra in the region mr 

Generally speaking, the existence of fixed points of RG-transformation does not guarantee 
that the critical asymptotics (^) do have place in the real system. As we have mentioned 
above, there would be another IR-divergences in the scaling functions JF^ and JF4 which 
are not handled by RG and, in principle, can modify the large-distance asymptotics close 
to the region mo <^ k. By the way, for mo/k 0, the critical dimensions are not sufficient 
to derive a conclusion on the long-range asymptotic behavior. To investigate the model 
(PPP with > in the region mo <^ A; in details we use a Short Distance Expansion 
method following f^. 

The short distance expansions of scaling functions JF^ and J-'a in (|47|) and (^81) provides 
us by an asymptotic relation of the form 

J^(l,mr) = l + ^Ci(r)m^[^'], (49) 

i 

where A[Oj] are the critical dimensions of all possible statistical momenta (the arbitrary 
products of fields and their derivatives averaged with respect to one point) of various 
quantities of the theory (0). In language of the modern critical phenomena theory 
such statistical momenta is called as composite operators in analogy with the well-known 
objects in quantum-field physics. 

It is obvious that the RG-predicted spectra (|4^) are still secure if for all operators 
A[Oi] > 0. The most important contributions into ( ^9]) for mr —>■ are those of the 
smallest A[0]. In the framework of e-expansion A[0] = dp + 0{e), where dp is the 
canonical dimension of O, therefore, if e is small, the canonical dimension do provides the 
major contribution to A[0]. That is why, in principle, to justify the scaling laws (^^, 
one can limit the checking of critical dimensions by the set of operators with minimal 
canonical dimension dg. At the leading order we consider the critical dimensions of the 
set of scalar quadratic operators with dp = Oi = 0^(x) and O2 = ^^(x). 

The certain critical dimensions are assigned to some linear combinations of the opera- 
tors Oi and O2 which still invariant in process of renormalization. The basis of renormal- 
ized composite operators are related to that one of non-renormalized operators through 
the renormalization matrix such that Fi = Zi^F^, 0. In principle, the calculation of 
matrix elements requires the analysis of diverging part of perturbation theory diagrams, 
however, the use of gauge symmetry consequences facilitates computations of the elements 

substantially. Since Oi is a gauge invariant operator, but O2 is not, the relevant Ward 
identities M prove the triangle structure for : Z22 = 1, Z21 = 0. One can say that 
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the non-invariant operator O2 does not contribute into scaling degrees of freedom of the 
gauge invariant theory, as well as it does not admix to the gauge invariant operator Oi in 
process of renormalization. 

Furthermore, we need not compute diagrams to determine the element Zu. Acting by 
the differential operation dm'^ onto the partition functional of renormalized theory (which 
is, obviously, finite with respect to the limit e 0), we obtain the following finite object 

(Z302)(x), (50) 

the finiteness of ( ^0] ) leads to the relation Zu = Z^^Zm = Z^. 

In principle, ( ^0|) give us enough information to compute the complete set of critical 
dimensions of the considered statistical momenta. Since Zj^ possesses a triangle structure, 
exactly the diagonal elements Zk^ give the relevant anomalous correction: 

e 

711 = 70 = g, 722 = 0. (51) 

Both linear combinations of Oi and O2 which have the definite scaling dimensions are 
also found unambiguously: 

Ci = Oi, A[Oi] = l + f, 
C2 = Oi + a02, A[02] = do = l. 

The only reason that we need to compute diagrams of perturbation theory for the element 
Z12 is to determine the value of a; this calculation gives a = 2 (this result is exact, i.e., it 
still valid for any value of e) 

For e > the most important contribution to ( ^Dj ) is provided by C2: ~ m, and as it is 
obvious, this contribution does not alter the scaling laws (^4|). The result (|5^) means that 
in the region mr the scaling laws (0) are still secure as the universal characteristics 
of the theory. 



7 Spectrum of Energy 

The proposed model provides a broad spectrum of practical results which can be com- 
pared with experimental data. We now compute the one-dimensional spectrum of energy 
of the fluid, E{r/lp), versus the dimensionless distance apart from the point of local per- 
turbation. As we have shown in the Introduction, the result on this spectrum, in principle, 
can be derived from phenomenology, (|T^. However, the question on justification of the 
phenomenological result ([T^ ) still remains, since the energy of perturbed fluid has two 
components, 

E{r/l,) ^ lmr){r/l,) + l{F'){r/l,), (53) 
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which would have dramatically different asymptotic behavior in the large-distance limit. 
By the way, would provide us an example of description of a conserved integral 
quantity in terms of scaling degrees of freedom. The spectrum (|53D is governed by the 
statistical momenta = ^{{d(f)Y){r /Ip) and E^ = ^{E'^){r /Ip) of the canonical dimen- 
sion (Ie = d. To determine the relevant critical indices we apply the trick which we have 
employed in the previous Section: Acting by the differential operations m^(9m2, gdg, and 
udu onto the partition functional of renormalized theory we obtain the following finite 
objects at the limit e — > 0: 



{Z,m'^'){r/lp), 

{[gdgZs] m202 + igOgZ,] {E^ + u^^A^) + i [4 + gd^Z,] g<P%r/lp), 
{[uduZ^] m'ct)^ + [ud^Z^] E^ + I [uduZ^] g<P^ + \ [ud^Z^] Ea+ 
+ [2Z^-uduZ2] u'ct?A^){r/lp). 



(54) 



Any linear combination of ( |54D is again finite at the limit e — 0, so that by means of 
simple arithmetic operations one can derive from ( |5^ ) the statements on finiteness for 
various linear combinations of statistical momenta including E^ and Ea- 

Taking into account that the scaling degrees of freedom relevant to scalar wave modes 
coupling are vanished in the large- distance limit {i.e., assuming (7 = in (|5^)), we obtain 
from (|5^ ) the combination containing E^ and E^'- 

{X,E^ + ^X2EA){r/lp), (55) 

where we have introduced Xi = [uduZi] and X2 = | [uduZ2]. In the framework of 
perturbation theory each of the coefficients Xi, as well as the momenta E^ and Ea have 
poles in e, and consequently each term in (^5|) separately is divergent if e — > 0. The 
meaning of (|55|) is that all the poles in e are subtracted out in such a way that the linear 
combination in ( ^5] ) has a definite limit for e — > 0. 
The linear form (|55| ) has two eigenvectors, 

Vi = E^ + X^'X2Ea, V2 = Ea + ^X^'XiE^, (56) 

such that the poles in e of Vi are subtracted out by those of Xi, and the poles in e of V2 
are eliminated completely by those in X2. By the way, = — -D^ In (uduZi) = — e/3 and 
IV2 = -D^ In iuduZ2) = -2e. 

For the real value er = 1/2 the major contribution to the energy spectrum for the 
large-distance asymptotics, r//p ^ 1, is provided by the combination from (|56| ) of minimal 
critical dimension, i.e., V2 with 



A[V2] = d + ^v2=d-l. 



(57) 
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Note, that A[Vi] = d + jvi = d — 1/6. Furthermore, for e,. = 1/2 X2 '^Xi = 2, and we 
have an exphcit form for Vi = + E^. The last step of the computation is to perform 



a Fourier transformation of the momentum asymptotics with the index (|57D into the real 
space {rlp^}. Bringing it about, we, finally, obtain a decaying profile (see Fig. 5), 

E{r/l,) ~ (^^^ (58) 
which meets the phenomenological result presented in the Introduction. 



8 Conclusion 

The final conclusion is that in contrast with either the statistical theory of waves (Za- 
kharov) or the theory of fully developed turbulence (Kolmogorov) the problem of coupling 
water to perturbation cannot be solved from phenomenology in principle. The matter is 
that the relevant physical system contains too many redundant degrees of freedom. That 
is why to fix the statistically stable behavior in the system one needs to add some extra 
assumptions on the character of perturbation as well as on the character of wave modes 
coupling. 

The problem of formulation of physically relevant hypotheses on wave modes coupling 
mechanism can be successfully solved by the use of various quantum field theory techniques 
and RG method in particular. This approach allows, first, to integrate over the redundant 
physical degrees of freedom, and, second, to investigate the asymptotic properties of 
physical systems by means of analysis of their scaling degrees of freedom. 

As a result we have formulated the effective action functional which allows to un- 
derstand the water response for a local external perturbation as a critical phenomenon. 
The critical system of water coupled to perturbation is subject to a "phase transition" 
depending on the certain physical properties of perturbation pulse. 

The results on the asymptotic behavior derived from the RG-analysis meet those of 
partial results which can be derived from phenomenology. 
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Table 1. The fixed points of the model 


N 


Coordinates {g'*,u'*} 


The stability region 


1 


{0,0} 


e < 


2 




unstable 


3 


{o,-^^} 


< e < 3 


4 


{26/3 + (2e)2|f,0} 


unstable 


5 




unstable 


6 


{-:^(i + yi-t(^/3-\/^),-y^} 


unstable 


7 


{-4(l-i/l-^(^/3 + y^)>V^} 


< e < 0.01 


8 


{-^(l-yi-^(^/3-V^)),-\/^^} 


3 < e < 3.35 
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CAPTIONS FOR FIGURES 



FIGURE 1. 



The potential energy U{f) versus /.a.) For the case of > the large-distance 
expectation value of / is trivial; b.) In the asymmetric case < the classical minimum 
is degenerate. There are two (in principle, an infinite number of) possible expectation 



values for /, ±y Starting from a given minimum, it is possible to describe all other 
minima by acting on the / with the U{1) symmetry group. 
FIGURE 2. 

Symmetrical and asymmetric phases of the model with no coupling to vector potential 
A; a.) > 0. The net large-scale momentum of the fluid outside the eddy which is 
formed around the perturbation region is equal to zero; b.) < 0. The pressure pulse 
rises the net fluid current into outside. P ^ 0. 

FIGURE 3. 

The diagram series of functions in the /^-model in the critical phenomena theory. 
FIGURE 4. 

The one-dimensional spectra for the pair correlation functions of the theory versus the 
dimensionless wave number. 
FIGURE 5. 

The one dimensional energy spectrum E{r/lp) versus the dimensionless distance apart 
from the perturbation point. 




